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ABSTRACT 


An automodel equation for the solar wind 
velocity is obtained in the self-consistent 
model. The solution of the convection- 
diffusion equation is obtained for the den- 
sity of galactic cosmic rays at a definite 
dependence of the diffusion coefficient and 
solar wind velocity on the rigidity and 
distance. 


In studying the propagation of galactic cosmic rays 
(GCR) in the interplanetary space, properties of the solar 
wind are considered, independent of the CR intensity. 

Though it should be noted that the influence of GCR on the 
supersonic solar wind is comparable with that of other 
factors (galactic magnetic field etc.). 

Then, while investigating the GCR modulations by the 
solar wind, one should take into account the fact that with 
the variation of the GCR intensity in time and space, the 
influence of GCR on the solar wind varies, which in its 
turn results in the variation of the solar wind velocity. 
That is, the self-consistent problem should be solved. 

In the spherically symmetrical case nonstationary 
equations of hydrodynamics for the solar wind with account 
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of the GCR influence have the form 
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where U , js are the velocity and the density of the solar 
wind, c a is the sound velocity, ft is the polytropic in- 
dex, Z is the heliocentric distance, and P is the GCR 
pressure equal to P - Ja/ (Z, R, t) pud. R (here A/ is the 
concentration of particles with the rigidity R ; p ; if are 
the momentum and the velocity of GCR particles, respective- 
ly). 

The variation of the spectrum N is determined by the 
character of GCR modulation in the interplanetary space, 
i.e. it is defined in the end by the velocity of solar wind 
and diffusion coefficient via nonstationary equation 

of diffusion with account of convection and GCR energy va- 
riation /1, 2/ 

If -v(*VA/)+(u v)N- (v u) -f ff- <2) 

Let us introduce as in /3/ new dimensionless dependent 

7 

variables 3 t and y by means of relations XL = — i 

c 2- * y -t*- and dimensionless independent variables Fni and 
En. Z » then the system (1) will have the form 

2C-oc + x £ +ac+-^-^'y+2y+y^£ tl pJ , J-f'yp' = 0 

(tnp)+ x ( £n.p) "+ x '+ 3-X =* O ^ 

y + xy' + 2 (x- 1 ) y+ 1 )y + 3?c ) = O 

where, in calculating the last term in the first equation, 
it is approximately assumed (ft- 1) ~ \ . 
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Let us assume that A and y are the function of 
one independent variable 2= i zt~ a ' , where a -const 
Let us require that Z and t do not enter explicitly 
in (2), then the density should have the form 
Since 3 /d(£ n t)-= -cl, )) 3 /a£ a the system 

(3) will take the form 

x , ("x-a 1 )+x 2 ~ac + -^(2y+-y' + y(Ea^) ^ + = q 

(En. $)(x- aj + a a + x' + 3x - O ^ 

V /y‘ ( < X-a f )+(&- 1)xt[a(vi) + (?fH)]x- 2 = o 

where the dash means a /a ZnZ , 

Eliminating d(0 a from the first two equations of 

the system, one may rewrite (4) as 

dEny 

cUn.2 = (*i-x) dx '(M) _ 

doc (3x-i)x-a 


(a 1 -cc) £ - 

oj f Q-a + + p ^ _ ocjj. x(l-x)(a 1 -ac) 

L ft 

where it is assumed F- dp /d^2 

As it is seen, the problem of integration of the ini- 
tial set of equations in partial derivatives of the first 
order is reduced to the integration of a first order ordi- 
nary differential equation of the form 


dy fi (x,-y) 

doc ' F a(x,X)) 

and to two quadratures. 

To define N one should solve eq.(2) which in the 
spherically symmetric case in dimensionless variables has 
the form 
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here S ) ^ Z A 0 , ^= r /r o , 'fc s Vt , a = T /? 0 , & - V?, 
where 1 a is the dimension of the modulating volume, R 0 
is the minimum rigidity of particles, t is the time 
which the solar wind spends to reach the boundary of the 
modulating volume. 

In order to obtain the analytical solution of eq.(6), 
one should choose the dependence of the diffusion coeffici- 
ent and the solar wind velocity in the form 

X = ot 0 / (, ~ 0() $ ZP ; u = u o/ £ "~^ 

where o( and are the arbitrary constants. 

Let us search for a solution with boundary conditions 
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Passing to the variable ^ ^ j=>°* and using the separa- 
tion of the variables sS (£ ; t ( 2 ) we have from 

( 6 ) 
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where K 2, is the separation constant. 
The solution of eq.(7) is 


«z) (8) 

where 3^ and y v are the Bessel functions of type 1 and 2, 
and "0 = f ~ —— . After satisfying the boundary conditions we 

2.<x. 

have 

Cy n+ ,(j«n)] a 
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